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We report a positive muon spin relaxation and rotation (pSR) study of the quantum spin ice 
materials Yb 2 Ti 207 and Yb 2 Sn 207 focusing on the low field response. In agreement with earlier 
reports, data recorded in small longitudinal fields evidence anomalously slow spin dynamics in 
the microsecond range below the temperature Tc at which the specific heat displays an intense 
peak, namely Tc = 0.24 K and 0.15 K, respectively, for the two systems. We found that slow 
dynamics extends above Tc up to at least 0.7 K for both compounds. The conventional dynamical 
Gaussian Kubo-Toyabe model describes the /rSR spectra recorded above Tc. At lower temperatures 
a published analytical extension of the Gaussian Kubo-Toyabe model provides a good description, 
consistent with the existence of short-range magnetic correlations. While the physical response of the 
two systems is qualitatively the same, Yb 2 Ti 207 exhibits a much larger local magnetic susceptibility 
than Yb 2 Sn 207 below Tc. Gonsidering previously reported ac susceptibility, neutron scattering and 
pSK results, we suggest the existence of anomalously slow spin dynamics to be a common physical 
property of pyrochlore magnetic materials. The possibility of molecular spin substructures to be 
associated to the slow dynamics and therefore the short-range correlations is mentioned. The slow 
spin dynamics observed under field does not exclude the presence of much faster dynamics detected 
in extremely low or zero field. 

PACS numbers: 75.40.-s, 75.10.Jm, 75.10.Kt, 76.75.-l-i 


I. INTRODUCTION 

The pyrochlore insulator compounds R2M2O'/, where 
i? is a rare earth ion and M a non magnetic element, have 
attracted considerable attention. The most remarkable 
discovery has been the spin-ice ground state of Ho2Ti2 07 
in 1997 A This ground state is now known for the emer¬ 
gence of a lattice-based Coulomb phase characterized 
by dipolar correlations and therefore pinch-point scat¬ 
tering patterns^ Its excitations are emergent magnetic 
monopoles.A^— 

This first success has led to systematic studies of other 
members of the R 2 M 20 g familyi^ However, apparently 
conflicting results have generated confusion. Only re¬ 
cently clarifying experimental data and thorough analy¬ 
ses have appeared. For example, low temperature mag¬ 
netization measurements for Tb2Ti207 barely support 
the existence of a magnetization plateau expected for a 
spin-ice like compound, at least down to 0.02 A 

first physical explanation of the persistent spin dynam¬ 
ics at low temperature as observed by the positive muon 
spectroscopy (/xSR) spin-lattice relaxation rate has been 
proposedi^ The dynamics would arise from low energy 
unidimensional excitations. They would be supported 
by spin loops originally discussed in Refs. [iMl 

Yb2Ti207 stands alone. The originally proposed low 
temperature spin liquid state, derived from powder sam¬ 
ple measurementshas been challenged ever since. 


with apparently no agreement reached on its basic phys¬ 
ical properties, such as the absence or existence of 
long-range magnetic correlations at low temperature as 
probed by neutron diffractionJAAAAI Remarkably, these 
correlations have recently been firmly established for a 
powder sample of the sibling compound Yb2Sn207A^A£ 

One of the most mysterious physical property of mag¬ 
netic pyrochlore compounds is the existence of anoma¬ 
lously slow spin dynamics. Here, a thorough study of 
the local magnetic fields and spin dynamics of Yb2Ti207 
and Yb2Sn2 07 as probed by /rSR spectroscopy under 
field is reported. Anomalously slow fluctuations of the 
Yb®“'" magnetic moments are detected below and above 
the temperature Tc at which specific heat displays a pro¬ 
nounced anomaly. Their effect on the measured spectra 
is fully described. Short-range correlations are detected 
below Tc for both compounds. Yb2Ti207 is found to be 
characterized by a much larger local magnetic suscepti¬ 
bility below Tc than Yb2Sn207. 

The organization of this paper is as follows. Sec¬ 
tion |TT] introduces the experimental method and explains 
the analysis of pSR spectra. In Secs. [Till and |TV] the 
Yb2Ti207 and Yb2Sn207 experimental spectra and their 
analysis are respectively presented. In the following sec¬ 
tion (Sec. El) a discussion of the present data and previ¬ 
ously reported extremely low-field results are provided. 
Conclusions are gathered in Sec. ED 
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II. EXPERIMENTAL AND DATA ANALYSIS 


The measurements were performed on powder 
samples already used for different experimental 
investigationsBoth of them are charac¬ 
terized by a well marked change in the longitudinal 
/iSR signal at Tc- For Yb^ Ti^O? this is in agreement 
with the results of Ref. [U but in contrast to Ref. [U 
for which a clear /iSR signature of only appears in 
the frequency shift. The magnetic transitions at Tc are 
first order with Tc ~ 0.24 and 0.15 K for Yb2Ti207 and 
Yb2Sn207, respectively. This was first shown in Refs, 
and [3 respectively. Later on, the first order nature of 
the transition for Yb2Ti207 was further characterized 
by detailed magnetization measurements.— Hence, the 
magnetic history of a compound is expected to influence 
results of measurements in its ordered state. When 
presenting our spectra for T < Tc we shall explicitly 
describe the experimental conditions under which they 
were recorded. Surprisingly enough, for both compounds 
no spontaneous /rSR magnetic field has been observed 
below Tc, even for a powder sample of Yb2Sn207 for 
which magnetic Bragg reflections have been detected as 

I 



mentioned in the introduction. 

The /rSR measurements were performed at the MuSR 
spectrometer of the ISIS facility (Rutherford Appleton 
Laboratory) and the LTF spectrometer of the Swiss 
Muon Source (S^S, Paul Scherrer Institute). 

The measurements were mostly done with the longi¬ 
tudinal geometry for which the external field Bext is set 
along the initial muon beam polarization In fact, 

it occurs that Bext and are not strictly parallel due 
to the presence of a separator in the beamlines. This is 
particularly true at LTF 3 .^ Hence, while the longitudinal 
polarization function derives from the counts detected 
in the detectors placed parallel and antiparallel to Bgxt, 
an oscillating signal is observed in the detectors set per¬ 
pendicular to Bext- Although its amplitude is small it 
contains valuable information which we have used in this 
study. The results obtained in this way were confirmed 
by recording a few spectra in the usual transverse geom¬ 
etry, i.e. with oriented perpendicular to Bgxti^ 

We now present the framework for the data analysis. 
The polarization function probing a static magnetic field 
distribution with field component distributions Dc{B°‘) 
identical for the three Cartesian directions a = A, Y, or 
Y, is derived from the Larmor equation. In the longitu¬ 
dinal geometry, it is expressed as 


Dc{B^)Dc{B^)Dc dB^AB^dB^, ( 1 ) 


with 52 = (sZ)2^ where 

7^ is the muon gyromagnetic ratio (7^ = 2 tt x 135.53 x 
10 ® rads“^T“^), and Hiong is the longitudinal field along 
the Y-axis which may not be necessarily equal to Text- 
The parameter Z describes a possible uniaxial anisotropy 
in the Z direction. Unless Z is specified, we shall assume 
Z = 1 , i.e. no anisotropy. 

In fact the field at the muon sites has often a dynamical 
character. The muon polarization function Pz{t) for a 
magnetic field fluctuating at rate Vc is computed using 
the strong collision model— > 2 ^ It satisfies the following 
integral equation; 


Pz(t) =P|*‘^*(t)exp(-zzct)-|- 


; f Pz{t - t')Pf’^^{t') exp{-iyct')dt'. 
Jo 


( 2 ) 


To proceed further an expression for Dc{B°') is needed. 
A conventional choice is a Gaussian function. Under¬ 
standing Pz{t) measured for Yb2Ti207 below Tc required 
the following simple extension of the Gaussian shape)^ 


P>c{x6) = ;^exp 


- ^imxf - ^imxf 


■ ( 3 ) 


Here, S sets the magnetic field scale and N is the nor¬ 
malization constant. In zero field and when % = 774 = 0 , 


i.e. for a Gaussian field component distribution, T|*’^*(t) 
is given by the famous Kubo-Toyabe formula— Under 
field and still in the Gaussian limit and with Z = 1 , 
the three-dimensional integral in Eq. [T] reduces to a one¬ 
dimensional integral— In zero field only symmetrized 
distributions are probed— In principle, in longitudinal 
field, /rSR can resolve asymmetric distributions. How¬ 
ever, concerning the data presented in this report, no 
clear deviation from symmetric distributions was found 
in a first analysis— Therefore, the distributions enter¬ 
ing Eq. [T] have always been symmetrized distributions 

£)sym(^0!) ag 

^sym(^a) ^ 1 [Bc{B^) + Pc(-P“)] - (4) 

Once T)®y™(P“) is determined, the standard deviation 
Alf of the distribution can be computed numerically ac¬ 
cording to 

/ OO 

b^Diy^ib)db, ( 5 ) 

-OO 

since the average of T)®^™(6) is 0 . If 773 = 774 = 0 , Alf = 
( 5 . 

Equation [ 3 ] is an ansatz which assumes the field com¬ 
ponent distribution to deviate only slightly from the 
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Gaussian shape. It can account for short-range spin- 
correlation effectsj^ The presence of short-range corre¬ 
lations is already known to transform the Kubo-Toyabe 
function into a polarization function with a shallower 
minimumi^ It is natural to apply Eq. [ 3 ] for the descrip¬ 
tion of longitudinal field (LF) ^SR asymmetry spectra 
recorded in a paramagnetic state. It can even be used 
below Tc when a spontaneous field is not detected, e.g. 
because of its dynamical naturei^ 

LF spectra will be analyzed with the two-component 
model 


[(I - F^,)Pz{t) + , (6) 

where is the initial /xSR asymmetry, Fbg the fraction 
of muons stopped outside of the sample, for which the 
relaxation is negligible. 

Except for some spectra recorded in the ordered state 
of Yb2Ti207, asymmetry time spectra corresponding to 
the transverse field (TF) geometry will be described with 
a weighted sum of two Gaussian damped oscillations: 




iTF 


(I - F2)exp - 


^2a2 +2 


F2 exp - 


y-2 


2 

cos(7^R2t + (k) 


cosij^Bit + (j)) 


( 7 ) 
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FIG. 1. (color online). fj,SK asymmetry spectra recorded at 
LTF in the paramagnetic phase of a Yb 2 Ti 207 powder sample 
under Bext = 18.8 mT. (left) Temperature scan of the longi¬ 
tudinal asymmetry. Circles are experimental data while solid 
lines represent fits as explained in the main text. The spectra 
for consecutive temperatures are vertically shifted by 0.1 for 
better visualization. A logarithmic time scale is used, (right) 
Thermal dependence of fitting parameters extracted for spec¬ 
tra recorded in the longitudinal and transverse geometries. 


A. Paramagnetic state 


Here, (j) is an experimental phase parameter and 
the initial asymmetry in the transverse geometry. For 
the vast majority of measurements for which LF and TF 
spectra were simultaneously recorded at LTF, ~ 
0 . 06 . The relation (Aq^)^ -|- (AJ^)^ = 0 . 27 ^ holds since 
the total asymmetry is a characteristics of the spectrome¬ 
ter. The mean fields Hi, B2, and the standard deviations 
Atf,i and Atf,2 refer to the first and second compo¬ 
nents, respectively. We find Atf,2 Atf,i, F2 = Fbg 
within experimental uncertainties, and B2 — Hext- These 
results are consistent with the attribution of the second 
component to the background, i.e. to muons missing the 
sample. Its narrow signal can serve as a precise mea¬ 
surement of Hext- For the first component we resorted 
to a simple function rather than trying to take into ac¬ 
count a possible deviation from the Gaussian field compo¬ 
nent distribution.^^ This is justified since, while we shall 
find deviations from the conventional description for the 
longitudinal asymmetries at low temperature, the simple 
expression in Eq. [ 7 ] provides a proper description of the 
measured transverse-field asymmetries. 


III. EXPERIMENTAL RESULTS: YbaTiaOz 

We shall first consider spectra recorded in the param¬ 
agnetic phase, i.e. above Tc ~ 0.24 K. In a second step 
data taken below Tc will be presented. 


We first show in Fig. [T] asymmetry spectra recorded 
with a longitudinal field Hext = 18.8 mT at different 
temperatures T > Tc- Remarkably, the asymmetry does 
not decay monotonically, but displays a well defined lo¬ 
cal minimum, at about 0.25 p,s at low temperature. The 
structure is less pronounced as the sample is warmed, but 
is still discernible at 700 mK, pointing to a relatively slow 
spin dynamics in the paramagnetic statei^ The spectra 
were analyzed with Eqs. [THBl The conventional dynam¬ 
ical Gaussian Kubo-Toyabe model was found to provide 
a good description, i.e. ??3 = 74 = 0. 

Transverse geometry spectra were analyzed with Eq. [ 7 ] 
and an example of such analysis is shown in the left panel 
of Fig. [71 In the right panel the Fourier transform (FT) 
amplitude of the spectrum is displayed. Only two peaks 
are resolved, in contrast to Ref. l 2 J where three peaks — 
one broad and two narrow — are reported for a powder 
sample, both above and below Tc. This will be further 
discussed at the end of Sec. iniBi 

The thermal dependences of the parameters extracted 
from the LF and TF analyses are presented in the three 
panels on the right of Fig. |T] As already noticed, B2 is 
a good measure of Hext- One can see that both Bi and 
Hiong are closely related and are substantially smaller 
than B2 — Hext- These important results will be dis¬ 
cussed in Sec. IV Al On increasing temperature this reduc¬ 
tion is attenuated. The fluctuation rate of the spin dy¬ 
namics is found to be i/c — 2 /xs“^ confirming the anoma- 
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FIG. 2. (color online), (left) Example of a transverse-field 
asymmetry spectrum measured at 300 mK under Bext = 
18.8 mT for a powder sample of Yb 2 Ti 207 . The spectrum 
was recorded at LTF with the spectrometer in the longitudinal 
geometry. This explains the small value for the asymmetry. 
Black circles display the spectrum while the red curve results 
from a fit to the data with Eq. [7] as explained in the main 
text, (right) The corresponding Fourier transform amplitude 
of the asymmetry data and fitting curve. 


lously slow spin dynamics already inferred from visual 
examination of the spectra. We find Alf and Atf,i to 
be roughly temperature independent. In contrast to ex¬ 
pectation, these parameters are quite different. However, 
we note that taking Z = 0 . 8 , the fit of the spectra are of 
similar quality and within uncertainties we still get the 
same results for Ba and r'c, but now Alf and Atf,i are 
approximately equal. The anisotropy of the field com¬ 
ponent distribution that this result suggests, could be a 
field-induced effect. The important point is that Bi and 
Siong remains nearly equal. 

In order to confirm the unusual slow spin dynamics we 
recorded data for a set of longitudinal fields at 300 , 450 , 
550 , and 700 mK, i.e. above Tc- Examples of spectra and 
fitting curves are presented in Fig. [ 3 l Again the dynami¬ 
cal Gaussian Kubo-Toyabe model gives a nice description 
of the spectra. Field-induced oscillations are clearly ob¬ 
served in a broad field range, supporting the existence of 
slow fluctuations. The fit parameters are plotted in Fig.| 4 ] 
as a function of Hext for temperatures up to 550 mK. The 
oscillations amplitude being somewhat reduced when the 
temperature is increased, some of the fit parameters are 
indeed poorly defined at higher temperatures. The re¬ 
duction of the field probed by the muons is confirmed; 
see the plot of B^ versus iJext- The parameters Alf 
and Atf.i characterizing the standard deviation of the 
Gaussian field distribution in LF and TF geometries, re¬ 
spectively, increase roughly linearly with Hext- 


B. Ordered magnetic state 

In a previous work it was shown that understanding a 
spectrum recorded at 200 mK, i.e. below Tc, and at low 
field for Yb2Ti207 required an extension of the conven¬ 
tional Gaussian distribution modelj^ This was confirmed 
with a numerical method based on the maximum entropy 
principle and the reverse Monte Carlo algorithm)^ In 



time t (/ 2 s) time t (fis) 


FIG. 3. (color online). Field dependence of asymmetry spec¬ 
tra recorded in the paramagnetic phase of a Yb2Ti207 pow¬ 
der sample in the longitudinal geometry. The data are rep¬ 
resented by black circles and the corresponding fitting curves 
as solid lines. The fits are explained in the main text. The 
spectra were recorded at LTF for T = 300 mK (left) and at 
MuSR for T = 700 mK (right). The data for consecutive 
Bext values (indicated nearby each spectrum) are shifted by 
0.1 unit for better visualization. 


YbaTbOj 



FIG. 4. (color online). External magnetic field dependence of 
physical parameters extracted from the analysis of field scans 
performed in the paramagnetic phase of a Yb 2 Ti 2 07 powder 
sample. We first plot Biong and Bi (collectively denoted as 
Be), then Vc, and finally Alf and Atf,i (collectively denoted 
as Ac). Most of the data have been derived from LF spectra 
(empty symbols), with few results obtained with the TF (full 
symbols) geometry at 300 and 550 mK. The plotting sym¬ 
bols correspond to different temperatures (given in units of 
millikelvin) and experimental geometries as indicated in the 
insets. 
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FIG. 5. (color online). Data recorded for a powder sample of 
Yb 2 Ti 207 at T = 70 mK with the longitudinal geometry, (a) 
Measured asymmetries with increasing Bext indicated nearby 
each spectrum. The spectra are vertically shifted by 0.05 for 
better visualization. The solid lines result from analytical 
fits as explained in the main text. For comparison, the blue 
dashed line represents the dynamical Gaussian Kubo-Toyabe 
fit of the asymmetry measured in a field of 10 mT. (b, c, d) 
Empty circles show the field dependence of three parameters 
describing the asymmetry spectra: the field Biong, the fluc¬ 
tuation rate and the standard deviation Alf = 1.70 5 for 
773 = 0.7335 and 774 = 0.4558. Full symbols represent the field 
dependence of B\, B 2 , and Atf,i obtained from the TF spec¬ 
tra (not shown). The dotted line in panel (b) corresponds 
to the curve Biong = Bext. (e) Field component distribution 
in reduced field scale. The solid line results from the anal¬ 
ysis with the extension of the Gaussian model. The dashed 
line is for the Gaussian field distribution corresponding to the 
dynamical Gaussian Kubo-Toyabe fit in panel (a). 


Fig. EKa) we display spectra recorded for different longi¬ 
tudinal fields at 70 mK, i.e. below T^, at LTF. Spectra 
taken with the same experimental conditions at MuSR 
were published elsewhere^ Here we profit from the bet¬ 
ter time resolution available at LTF. The sample was 
zero-field cooled. The data were subsequently recorded 
for increasing values of the field. Weak oscillations are 
observed, with an amplitude which is again reduced at 
high field. The asymmetry spectra were analyzed with 
Eqs. [IHBl The separate analysis of the spectra led to 
parameters 773 and 774 independent of the field, within 
uncertainty. Therefore in a second step, the spectra at 
different fields were simultaneously fit with the same set 
of parameters 773 and 774. The fits shown in Fig. [SKa) are 
excellent. We found 773 = 0.73 ( 2 ) and 774 = 0.46 ( 2 ). This 
is consistent with the previous analysis of a spectrum 
recorded at 200 mK under 2 mT.^^ The two parameters 
773 and 774 are not negligible. This means that 
displayed in Fig. Ue) deviates substantially from the con¬ 
ventional Gaussian shape as the comparison shows. From 


7 “-I 

H 

g 0.08 

0.06 
3 

I 0.04 
H 

■ 0 1 2 3 

time t {fis) magnetic field (mT) 

FIG. 6 . (color online), (left) Example of a transverse-field 
asymmetry spectrum measured at 50 mK under Bext = 50 mT 
in field cooling condition for a powder sample of Yb 2 Ti 207 . 
The spectrum was recorded at LTF with the spectrometer in 
the transverse geometry. Black circles display the spectrum 
while the red curve results from a fit with a model similar to 
Eq.[7]but with three Gaussian oscillating components, (right) 
The corresponding Fourier transform amplitude of the asym¬ 
metry data and fitting curve. The broad signal stems from 
the sample while two narrow peaks at B^ ~ 46 mT and 
B^ ~ 49 mT originate from muons stopped outside the sam¬ 
ple. 
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the values of 773 and 774 we compute the field distribution 
standard deviation Alf = 1.70 ( 9 ) 5 . Note that the pa¬ 
rameters 773 and 774 are strongly correlated. The inferred 
D®^™(i?“) is characterized by relatively pronounced tails. 
The fitting parameters Biong, I'c, and Alf are presented 
in Figs. [ 5 Kb)-[ 5 Kd). We find T^long ^ 7 ?ext with Biong 
roughly proportional to Bext, as in the paramagnetic 
state but with a much smaller coefficient of proportion¬ 
ality. The fluctuation rate z/c ^ 1 is smaller than in 
the paramagnetic state but still appreciable. The stan¬ 
dard deviation Alf increases quite linearly with field, in 
agreement with the behavior in the paramagnetic state. 
Atf,i is slightly different from Alf at high field. As 
for the paramagnetic state, this difference becomes neg¬ 
ligible if a small field-induced anisotropy is introduced. 
This is a side effect which can be neglected, as already 
mentioned when the paramagnetic data were discussed. 

The strong reduction of Biong relative to Bext is re¬ 
markable. As expected, and supporting our data anal¬ 
ysis, Bi ~ Biong. We expand on the meaning of this 
reduction in Sec. IV Al 

In order to resolve the discrepancy between our re¬ 
sults and those reported in Ref. we performed addi¬ 
tional measurements at the LTF spectrometer in TF con¬ 
figuration. This allows us to simultaneously record LF 
and TF spectra with a large = 0 . 21 , and a smaller 
Aq^ = 0 . 17 . The experimental conditions, namely mea¬ 
surements in field cooling protocol with Bext = 50 mT 
and at T = 50 mK, are as in Ref. The result is shown 
in Fig. [6] Qualitatively our Fourier spectrum is con¬ 
sistent with that reported in Ref. [l^. Three peaks are 
visible. This is in contrast to the spectrum displayed in 
Fig. m for which only two peaks are present. Regarding 
this latter spectrum, Bext is smaller and the data concern 
the paramagnetic instead of the ordered state. While it 
is natural to assign the broad peak around 20 mT in Fig. 
[S]to the sample, we believe that the other two arise from 
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FIG. 7. (color online). Results obtained in the paramag¬ 
netic phase of a powder sample of Yb 2 Sn 207 under Rext = 
18.8 mT. (left) Examples of asymmetry spectra obtained at 
LTF in the longitudinal geometry. Circles are experimental 
data while solid lines represent fits as explained in the main 
text. The temperature is given nearby each of the spectra 
which are vertically shifted by 0.1 units for better visualiza¬ 
tion. (right) Parameters extracted from a ht to Eqs.[TH6]for 
LF data and Eq. [7| for TF data. 



FIG. 8. (color online). Examples of field scan measurements 
of asymmetry spectra performed in the paramagnetic phase of 
a powder sample of Yb 2 Sn 207 (black circles) and correspond¬ 
ing fitting curves obtained with Eqs. [T}{6l The spectra were 
recorded in the longitudinal geometry at the LTF spectrome¬ 
ter for T = 250 mK (left) and T = 550 mK (right). The Rext 
value is given nearby the spectra which are shifted by 0.1 unit 
relative to each other for better visualization. It is remarkable 
that oscillations due to Rext are clearly observed pointing to 
a low frequency spin dynamics above Tc in Yb 2 Sn 207 . 


muons stopped outside the sample. Indeed, the relative 
contributions of the peaks at 46 and 49 mT to the Fourier 
spectrum are 14(1)% and 12(1)%, respectively. Assign¬ 
ing only the peak at 49 mT to the background leads to an 
unreasonably small background value for the LTF spec¬ 
trometer. Most probably, the stray field and the geome¬ 
try of the sample and sample-holder are at the origin of 
the two peaks on the right of the field scale. The narrow 
peak should correspond to the part of the sample-holder 
located far from the sample while a slightly broader peak 
originates from the sample-holder part in the immediate 
neighborhood to the sample. 

In any case, the agreement demonstrated in Fig. [51(b) 
between the mean fields i?i, measured directly from TF 
spectra, and Siong, extracted from the analysis of LF 
spectra, is a strong support for the methodology of our 
TF and LF data analysis. 


IV. EXPERIMENTAL RESULTS: Yb2Sn207 

Here we present results of the analysis for Yb 2 Sn 2 07 
in the same way as was done for Yb 2 Ti 207 . We first 
focus on the paramagnetic state. 


A. Paramagnetic state 

In Fig. [7] we display results obtained in a LF of Hext = 


18.8 mT above Tc = 0.15 K. Also in this compound one 
can clearly see oscillations in the asymmetry spectra as¬ 
sociated with Hext, which directly point to a slow spin dy¬ 
namics. The dynamical Gaussian Kubo-Toyabe model, 
i.e. Eqs. (THU with 773 = 774 = 0 provides a good account of 
the spectra. The local mean field, i.e. B\c,ng, in this com¬ 
pound is again reduced. The strength of the reduction 
is weaker than in Yb 2 Ti 207 and it also decreases with 
increasing temperature. The difference between Bi and 
Siong is small and may not be really significant. The spin 
fluctuation rate i>c gradually increases with temperature 
and reaches a maximal value of 2.0 (2) /is“^ at 700 mK. 
As expected, the magnitudes of Atf.i and Alf are al¬ 
most equal. This is in contrast to Yb 2 Ti 207 . 

In order to clarify the origin of the slow spin dynam¬ 
ics in the paramagnetic state additional measurements at 
250, 420, and 550 mK were performed for different fields. 
In Fig.jSjwe show examples at 250 and 550 mK, together 
with the fitting curves depicted as solid lines. The oscil¬ 
lations due to Hext are present in the whole field range of 
the measurements. Again, a Gaussian field component 
distribution is sufficient for the description of the spec¬ 
tra. In Fig. 0 the corresponding fitting parameters are 
presented as a function of field. The field dependence of 
the three parameters measured for the stannate and the 
titanate in their paramagnetic states is rather similar. 
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FIG. 9. (color online). Same caption as for Fig. [d] Here 
the physical parameters concern the paramagnetic phase of 
YbaSnaOr. 


B. Ordered magnetic state 



Spectra were recorded in the magnetically ordered 
state at 14 mK. They were presented in Ref. m. The 
sample was first zero-field cooled and then Bgxt was ap¬ 
plied with the following sequence: 5, 3.4, 9.7, 19, 50, 
and 100 mT. The LF and TF asymmetry spectra are 
well described with Eqs. [THB] and Eq. [3 respectively. As 
for Yb 2 Ti 207 the parameters and 774 are appreciable 
and field independent within the statistical uncertainties. 
Therefore, we performed a global fit of the LF spectra 
with common 773 and 774 in order to determine the field 
distribution. The experimental LF time asymmetries to¬ 
gether with the fitting curves are presented in Fig. unia). 
The analysis results in 773 = 0.67 (4) and 774 = 0.40(3), 
i.e. only slightly different from the values measured for 
Yb 2 Ti 207 . From these values we get Apr = 1.38(9)5. 
In Figs. [T0l' b)- fl0r d) we show the field dependence of the 
other LF parameters, together with the TF parameters. 
A linear dependence of the field standard deviations is 
also present in the ordered state. The fluctuation rate Vc 
is slightly reduced and exhibits a weaker field dependence 
than in the paramagnetic state. The fields i?iong and Bi 
are equal within experimental uncertainties. They are re¬ 
duced relative to Bext, but not as strongly as in the case 
of Yb 2 Ti 207 ; see Fig. [5] In Fig. [rof e) we compare the 
results of the fits with the Gaussian model and with its 
extension. The latter solution provides a better descrip¬ 
tion. The field component distribution obtained from 
this analysis is shown in Fig. fTOT f) as a solid line. As 
for Yb 2 Ti 207 , a significant deviation from the Gaussian 
shape displayed as a dashed line is evident: pronounced 
tails are present. 


FIG. 10. (color online). Data recorded at the LTF spec¬ 
trometer for a powder sample of Yb 2 Sn 2 07 at T = 14 mK 
with the longitudinal geometry, (a) Field dependence of LF 
time spectra (circles) and fits (solid lines) as described in the 
main text. The spectra are vertically shifted by 0.07 for bet¬ 
ter visualization (b, c, d) Corresponding fitting parameters as 
a function of Rext. The field component standard deviation 
is Alf = 1.38 5 for 773 = 0.670 and 774 = 0.400. The Bext 
dependence of the parameters deduced from the TF spectra 
is also displayed, (e) LF time spectrum (circles) at 3.4 mT 
plotted with a linear time scale and fits with the extension 
of the Gaussian model in red solid line and the conventional 
Gaussian Kubo-Toyabe model in blue dashed line, (f) Field 
component distributions in normalized field units. The solid 
line shows the distribution obtained from a global fit using 
the extension of the Gaussian model, and the dashed line the 
Gaussian distribution corresponding to the Kubo-Toyabe fit 
in panel (e), also in blue dashed line. 


V. DISCUSSION 


We shall first consider the conclusions reached here 
from the analysis of spectra recorded under longitudi¬ 
nal fields. Then we shall discuss their compatibility with 
those previously derived from extremely low field mea¬ 
surements. 































A. Measurements under relatively large fields 

Let us first summarize the results of the analysis of 
the /iSR spectra for the two compounds discussed in this 
report, i.e. Yb2Ti207 and Yb2Sn207. 

The following features emerge. From the highest tem¬ 
perature investigated, i.e. 700 mK, down to far below 
the temperature at which a sharp peak in the specific 
heat is detected, the fluctuation rate of the Yb^+ mag¬ 
netic moments is equal to ~ l/is“^, independent of 
the field and temperature within a factor 2 . Hence, the 
spin dynamics is anomalously slow. The component field 
distribution measured under longitudinal field is Gaus¬ 
sian above and somewhat deviates from it below T^. 
Its standard deviation increases linearly with the field 
for the two compounds, both below and above Tc- When 
normalized to the width, the distribution is field indepen¬ 
dent in the probed held range. Finally, a reduction of the 
mean held at the muon site is measured for the two com¬ 
pounds. Within experimental uncertainties, it is similar 
in longitudinal and transverse held measurements. It is 
particularly strong for Yb2Ti207 at T < Tc. Remark¬ 
ably, it is approximately linear in held, both below and 
above Tc. The ratio (Rq, — Rext)/Rext can therefore be 
considered as a local magnetic susceptibility. 

As the magnetic phase transitions are hrst order, the 
procedure followed to record spectra below Tc has an 
inhuence on the results of the measurements. We recall 
that the asymmetry spectra analyzed in this report were 
taken after cooling the samples in zero held, except for 
the data of Fig. | 6 l 

We now discuss the meaning of our results. We hrst 
note that the existence of slow spin dynamics in mag¬ 
netic pyrochlore compounds is well documented. It was 
reported in 2002 for Yb2Ti207 from the analysis of /rSR 
and Mbssbauer measurementsSoon after, it was re¬ 
alized that a wide spectrum of huctuation rates exists 
in this compoundji^ Later on, slow spin dynamics was 
reported for the pyrochlore Tb2Ti207 using ac suscepti¬ 
bility measurements!^ Slow spin dynamics with a wide 
range of huctuation rates of the rare-earth magnetic mo¬ 
ments in geometrically frustrated magnetic compounds 
is in fact one of their characteristics. As an additional 
example, we mention the garnet compound Gd3Ga50i2 
with three-dimensional hyperkagome magnetic lattice 
for which neutron-scattering studies reveal a remarkable 
wide range of energies from milli- to picoelectronvolt.— 
The present study puts on hrm ground the existence of 
slow dynamics for Yb2Ti207 and Yb2Sn207 as seen by 
/iSR. 

As to the physical mechanism at the origin of the slow 
dynamics, we note that persistent spin dynamics at low 
temperature, i.e. below ~ 1 K, as observed by ^SR spin- 
lattice relaxation in pyrochlore compounds has recently 
been proposed to arise from unidimensional excitations 
that would be supported by spin loopsi^ These one¬ 
dimensional objects in three-dimensional crystal struc¬ 
tures such as the pyrochlore structure have been studied 


early ouj ^^d^ Their existence provides a natural expla¬ 
nation for slow dynamics since it means that spin loops, 
i.e. a large number of spins in contrast to a conventional 
single spin, are involved in the dynamics. This picture 
may extend above ~ 1 K, the interaction between spin 
loops being recognized. This is in qualitative agreement 
with the observation of slow spin dynamics deep into the 
paramagnetic state of Nd2Sn207i^ Interestingly, classi¬ 
cal simulations of spin dynamics for the antiferromag¬ 
netic Heisenberg model on the kagome lattice is consis¬ 
tent with the loop picture; see Ref. IH and references 
therein. 

The field component distribution extracted from the 
measurements below Tc deviates from the Gaussian shape 
for both compounds, pointing out the existence of short- 
range correlations. These correlations are absent in 
the ordered magnetic state of a conventional compound. 
They do exist for Er2Ti207i^ The increase of the field 
standard deviation with Text suggests that the held in¬ 
duces some disorder. The mechanism at play is not clear. 
It could be that the held induces a canting of the mag¬ 
netic moments. Surprisingly, we have not found any 
effect of magnetic short-range correlations on the held 
component distribution in the paramagnetic state. Yet 
such correlations, signalled by a broad hump in the spe- 
cihc heat centered at about 3 and 2 K for Yb2Ti207 and 
Yb2Sn207, respectively, do existi^i^iiiSiiL and have even 
been observed by neutron scatterin g^^d^ in a crystal of 
Yb2Ti207 whose specihc heat does nonetheless display 
no sharp peak at Tci^i The detection of these correla¬ 
tions by p,SR needs to be further investigated. 

Let us now discuss the origin for the reduction of the 
mean held at the muon site. We have the key result that 
the same reduction is measured for LF and TF asymme¬ 
tries. This suggests that in fact we are dealing with a 
frequency shift effect and not an effect of intermittency 
of the held at the muon sitei^ This interpretation is sup¬ 
ported by the results published for Tb2Ti2074^ In this 
system, below 10 K, the relative frequency shift varied 
between — 0.2 and — 0 . 7 . In Yb2Ti207 [Fig. [ 5 l[b)] it is 
approximately — 0.6 and for Yb2Sn207 [Fig. ITUI bl] it is 
near — 0 . 2 . In conclusion, there is no exotic physics in 
the measured held reduction. It just occurs that the fre¬ 
quency shift is negative. If the muon site localization 
were known, it would probably be possible to understand 
this feature. This has recently been done for MnSi 4 ^ The 
much larger reduction below Tc for Yb2Ti2 07 relative to 
Yb2Sn207 is just the signature of a much larger local 
magnetic susceptibility for the former compound. 


B. Compatibility with previously published 
extremely low field measurements 

The asymmetry and therefore the polariza¬ 

tion function Pz{t) in the paramagnetic state does not 
decay monotonically at intermediate Rext values for 
Yb2Ti207 and Yb2Sn207. On the other hand, at ex- 
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tremely low fields it doesJ^d^ Quantum perturbation 
theory describes Pz{t) in terms of correlation functions 
when the spin dynamics is sufficiently fast,— namely 
i^c/(7/iALF) ^ 2 — It has been shown that this quantum 
description and the one exposed in Sec .HIl give consistent 
results in this limit when the field component distribu¬ 
tion is Gaussian)2& as in the paramagnetic state of the 
two compounds of interest. Remarkably, at extremely 
low Rext the quantum description leads to a monotonic 
decay of Pz{t) which is exponential-like for P 1 — 

In our case, the shortest time we probed is t ~ O.I /xs. 
This means that i^c ^ 10 a lower bound larger than 
the value extracted from the analysis of measured 

under field which is « 1 /is“^. Hence, the fluctuation 
modes probed at extremely low field are fluctuating much 
faster than at intermediate field. 

This discussion gives results entirely consistent with 
previous inferences from extremely low Hgxt measure¬ 
ments near Tc in the paramagnetic state— From 
i 70 Yb Mossbauer measurements the fluctuation rate for 
Yb2Ti207 and Yb2Sn207 was found to be in the 10^° s“^ 
range. Identifying ^SR and Mossbauer fluctuation rates, 
taking into account the measured spin-lattice relaxation 
rate Xz characterizing the decay of Pz{t), and the con¬ 
ventional motional narrowing formula Xz = 27^Alp/7'c, 
Alf was computed to be consistent with expectation for 
the known Yb^’*' magnetic moments. 

At this juncture we note that the field standard de¬ 
viation related to the slow dynamics is more than one 
order of magnitude smaller than the one characterizing 
the faster dynamics. This means that only a small frac¬ 
tion of the Yb^+ magnetic modes is involved in the slow 
dynamics. 

Surprisingly, no spontaneous /xSR oscillation is ob¬ 
served below Tc for the stannate. Such a lack of oscilla¬ 
tion was also observed for Tb2Sn207 and interpreted as 
a signature of the dynamical nature of the ground state 
of the compound,later on confirmed by neutron-spin 
echo measurements— This means that Vc > 27 ^i 3 fluc, 
where Vc is the fluctuation rate of the fluctuating spon¬ 
taneous field Bfluc at the muon site— From the relax¬ 
ation rate Xz — 2 . 3 us“^ measured at low temperature 
in Tb2Sn207 (Ref. l 34 h and from the fluctuation rate Pc 
= 5 X 10 ^*^ s“^ measured with the neutron spin echo 
technique, we deduce Hfluc = 0.4 T using the formula 
= 7 ^^iuc/^c- Scaling Rfluc to the ordered mag¬ 
netic moments — 5 . 4 /xb in the Tb2Sn207 and approxi¬ 
mately l.I /iB in the Yb pyrochlore o^^d^ — we compute 
Bfiuc — ( 1 . 1 / 5 . 4 ) X 0.4 = 0.08 T for the two ytterbium 
compounds. Therefore Pc > 7 x 10 ^ s“^. 

Hence, the analysis of the fjSR data leads to two 
magnetic fluctuation rates in the paramagnetic state: 
Ri 10® s“^ and ~ 10^® s“^ close to Tc. In the ordered 
state two rates are also detected: again ss 10® s“^ and a 
rate larger than 7 x 10 ^ s“^. Finding more than one fluc¬ 
tuation mode is not surprising: we refer to the ordered 


state of Tb2Sn2 07-^‘^^-^ 

VI. CONCLUSIONS 

In this paper we have presented a complete /xSR study 
of the spin dynamics of Yb2Ti2 07 and Yb2Ti2 07 as seen 
in a modest applied field. Extremely slow fluctuation 
modes with a scale of approximately 10® s“^ were found 
both below and above the transition temperature. On 
the other hand, a much faster fluctuation rate, namely 
in the 10^® s“^ range near Tc, was recognized from ex¬ 
tremely low field paramagnetic measurements. In addi¬ 
tion, a scale larger than 7 x 10 ^ s“^ does exist in the 
ordered state. A similar spectrum of fluctuation rates is 
present in the two compounds. Such a broad and complex 
spectrum of fluctuations is unusual in a magnetic system. 
To our knowledge no theory is available for comparison. 

The field component distribution deviates from the 
Gaussian shape for the two compounds below Tc, point¬ 
ing out the existence of short-range magnetic correla¬ 
tions. This is an unconventional feature which suggests 
the formation of molecular spin substructures, such as 
spin loops. They probably also exist in the paramag¬ 
netic state, although a signature of them has only been 
found from specific heat. Finally our study reveals that 
the magnetic susceptibility of Yb2Ti207 is much stronger 
than that of Yb2Sn207 below Tc. 

Still considering the ytterbium pyrochlore family, we 
note that ferromagnetic superexchange interactions drive 
the magnetic transition observed at Tc for Yb2Ti207 and 
Yb2Sn207. These compounds occur to be spin-ice like. 
It would be of interest to investigate with /xSR techniques 
Yb2Ge2 07 for which antiferromagnetic superexchange 
interactions are responsible for the magnetic transition— 
In particular, possible differences in the spin dynamics 
would be of interest. In this respect, we note that the 
relaxation for the antiferromagnet Er2Ti207 is strikingly 
different from observed in this report] ®®:®® In particular, 
the field responses below Tc are quite different. 
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